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ABSTRACT 


111 


The United States National Aeronautics and Space Administration in 
cooperation with 12 scientific groups launched 26 sounding rockets from the 
equatorial launch site at Punta Chilca, Peru during the period May 20, 1975 
to June 10, 1975. The purpose of this cooperative effort is to study in 
detail the structure and dynamics of the equatorial atmosphere between 20 
and 160 km altitude. Three of the 26 rocket payloads were buil - ' by the 
Aeronomy Laboratory of the Department of Electrical Engineering of the 
University of Illinois at Urbana-Champaign . These payloads were designed 
to investigate the anomalous properties of the equatorial ionosphere. The 
objectives of the UI payloads included 

1. A measure of electron temperature, 

2. A measure of electron concentration in the lower ionosphere, 

3 . An observation of the fine structure in the profile of electron 
concentration, and 

4. An examination of the role of energetic electrons as a nighttime 
source of ionization. 

This report discusses in detail a technique for measuring electron 
concentrations in the lower portion of the ionosphere above Punta Chilca. 
The technique combines a radio -propagation experiment for measuring Faraday 
rotation and a dc/ Langmuir probe experiment for measuring electron current . 
The results obtained from the analysis of radio and probe data from Nike 
Apache 14.532, which was launched at 20:26 UT on May 28, 1975, at a solar 
zenith angle of 60°, are presented in Chapter 5. A comparison of the 
profiles of electron concentration, N, from Nike Apache 14.532 and Kane 
[1974] indicates that the value of the maximum ionization in the D region, 
i!7 D, under quiet conditions is proportional to the square of the cosine of 
the solar zenith angle. 
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1. INTRODUCTION 

Twenty-six sounding rockets were launched during the period May 20 to 
June 10, 1975 by the United States National Aeronautics and Space Adminis- 
tration in cooperation with 12 scientific groups to study in detail the 
structure and dynamics of the equatorial atmosphere between 20 and 160 km 
altitude. The rockets were launched from Punta Chilca, Peru, at the geo- 
magnetic equator. 

The’ Aeronomy Laboratory of the Department of Electrical Engineering 
of the University of Illinois at Urbana-Champaign built three of the 26 
rocket payloads. The payloads were designed to investigate the anomalous 
properties of the equatorial ionosphere. More specifically, the objectives 
of the UI payloads included 

1. A comparison of electron temperatures measured by three techniques 
(a) The University of Illinois Langmuir probe 

Cb) An RF resonance probe supplied by Dr. K. Hirao and Dr. K. 

Oyama of the University of Tokyo 
Cc) The incoherent-scatter radar at Jicamarca, Peru 

2. A measure of the electron concentration in the lower ionosphere ; 

3. An observation of the fine structure in the profile of electron 
concentration to examine the role of vertical transport in layer- 
ing metallic ions; and 

4. An examination of the role of energetic electrons as a nighttime 
source of ionization. 

This thesis discusses the technique used for measuring the electron 
concentrations in the lower portion of the ionosphere above Punta Chilca. 

The technique combines a radio -propagation experiment for measuring Faraday 
rotation and a Langmuir dc -probe experiment for measuring electron current . 
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Although the probe cannot measure electron concentration directly, its high- 
resolution measurements of changes in electron concentration are calibrated 
by the radio -propagation experiment which accurately measures electron con- 
centration but with poor altitude resolution. 

A brief discussion of the Appleton-Hartree and Sen-Wyller theories is 
presented in Chapter 2. Chapter 3 describes the operating principles of the 
radio -propagation and dc -probe experiments. A method of processing the 
radio-propagation data to obtain Faraday rotation (differential phase) rates 

and a method, based on the generalized magnetoionic theory for analyzing the | 

Faraday rotation rates is presented in Chapter 4. These methods are applied j 

to radio and probe data from Nike Apache 14.532, which was launched at j 

20:26 UT on May 28, 1975, at a solar zenith angle of 60°. The data from the j 

' ' ■ . • • • . i 

two nighttime shots, Nike Apaches 14.524 and 14.525, are still being analyzed * • 

and will not be presented here. 

- ; . ' • ; . • ' ^ 7 -j 

■■ ' J j 

. ' . ... . \ 
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2 . MAGNETOIONIC THEORY 


The equations of magnetoionic theory describe the propagation of 
radio waves through ionized media, upon which is superposed an external 
magnetic field. The ionosphere and magnetic field of the earth constitute 
a magnetoionic medium. In upper regions of the ionosphere, collisions of 
electrons are negligible, and magnetoionic theory reduces to its simpler 
form, the Appleton-Hartree equations [Appleton, 1932; Eartvee, 1931] which 
are useful for discussion of wave polarization. In the lower region of the 
ionosphere, collisions are important. They result in the conversion of 
electromagnetic energy from the radio waves to heat. In this region a more 
generalized theory is necessary. Generalized magnetoionic theory is ex- 
pressed by the Sen-Wyller equations [Sen and Wyller, 1960] which must be used 
in discussions of radio-wave absorption. The Sen-Wyller equations are used 
in Chapter 4 for the analysis of radio -propagation data. 

2.1 Appleton-Eartvee Theory 

An important property of wave propagation in an anisotropic medium, 
such as the earth's ionosphere, is the occurrence of normal modes or 
characteristic waves. A characteristic wave is a wave which propagates in 
a given uniform medium without changing its wave polarization. In an aniso- 
tropic medium there exist two characteristic waves, each having a distinct 
polarization and refractive index. The theory of Appleton-Hartree describes 
the polarizations and the refractive indices of these modes in the ionosphere 
under the following assumptions: electron collisions with neutrals are in- 

dependent of electron energies, the medium is electrically neutral with a 
uniform charge distribution, the magnetic field is uniform throughout the 
medium and the wave frequency is much larger than all ionic gyrofrequencies . 
When the wave frequency is much larger than all the ionic gyrofrequencies, 


I 
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induced ionic motions are negligible because of the relatively large mass 
of the ions. Thus, the ions form a stationary neutralizing background. 
Maxwell's curl equations for plane waves of the form 
E = E n exp U (ut-tz)] 


are 


D = {-n/o)z x 
H = (e 0 /y 0 )' 5 nz x E 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where 


k = wave number (rad m - -*-) 
w = wave frequency (Hz) 
n = index of refraction 

0 = velocity of light (m/s) 

= permittivity of free space (F/m) 
\i 0 = permeability of free space 
t = magnetic field intensity 

1 = electric flux density {C/m) 

E = electric field intensity 


and 


E 


E^x + E u y 


o J ixr ' “y* 

The wave polarization p is defined [Budden, 1961] as 


P - 


H E P 

_£ = JL = JL 

H E P 

y x x 


(2.4) 


(2.5) 


where 


P r and P are the polarization densities {C/m ) in the x and y directions. 

y 

The indices of refraction for fields in the form of equation (2.4) are given 
by 


i* 1 

n = 1 + 




e c>Ex 


- 1 + -**■- 
e J? 


O y u/e7u 


( 2 . 6 ) 


o»o 


1 


J 
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The electric field intensity is related to the polarization density by: 
zE = T 1 * P (2.7) 

where the inverse of the susceptibility tensor is 


7-1 -1 

* = T 


l -tf L ji f ' 




0 

1 


for regions in which collisions may be neglected, and 


to 


X = 




co 


Ne^ 

z Q m 


( 2 . 8 ) 


„ -eB 

Y t = — — cos 
L m 


,, -eB . , 

*T " — Sin * 


co^ = plasma frequency (rad/s) 
e = magnitude of electronic charge (C) 
m = electron mass (kg) 

_3 

iV 5 electron concentration (m ) 

B = magnitude of the magnetic flux density, ^ (Wb/m^) 

(j> = angle between the wave normal and B 
Combining equation (2.7) with Maxwell's equations (2.2) and (2.3) results 
in the following quadratic equation in p 


p' - J-jzjr p + 1 0 


(2.9) 


Applying the quadratic formula to (2.9) yields the following roots : 

2 . 

' L 1 

P = 3 


V /2 A . 


1-X \ 1+J / 

The solution for the square of the refractive index is 


n = 1 - 


1+JI L P 


(2.10) 


(2.11) 
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Equations (2.10) and (2.11) are known as the Appleton-Hartree equations 
(collisionless case) . 

For rocket flights from Punta Chilca, Peru (Geog, Lat. 12° 30' S, 

Geog. Long. 76° 48' W) the direction of the propagated radio wave is ap- 
proximately transverse to the earth's magnetic field. In this case, the 
quasi-transverse approximation of the Appleton-Hartree equations is valid. 

l//4Y L 2 » |(1-X) 2 | (2.12) 


Applying this approximation to equations (2.10) and (2.11) yields the 
following results: 



(2.13a) 

(2.13b) 


p x ~ ” 

n ^ = 1 

x 


X 


(2.14a) 

(2.14b) 


where the subscripts o and x refer to the ordinary and extraordinary waves, 
respectively. Equation (2.13a) shows that the ordinary wave is linearly 
polarized with its electric vector parallel to the earth's magnetic field, 
while equation (2.14a) shows that the extraordinary wave is linearly polar- 
ized with its electric vector perpendicular to the earth's magnetic field. 
2.2 The Sen- Wylie? Theory 

In the region of the lower ionosphere between about 50 and 90 km, the 
collisions of electrons with neutral molecules appreciably influence the 
refractive and absorption indices of radio waves which are suitable for in- 
vestigating this region. The proportionality of electron collision frequency 
and electron kinetic energy established by Phelps and Pack [1959] is an 
essential part of the equations for the numerical analysis of experimental 
data. The generalized magnetoionic theory which takes into account the 
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proportionality of electron collision frequency and electron kinetic energy 
was developed by Sen and Wyller [1960] . 

Sen and Wyller expresed the refractive index and wave polarization 
as functions of the elements of the permittivity tensor, e. 


G4+Bsin 2 ij>±(B sec %-C 2 cos 2 <t>) ) 
n 2 

D + E sin <j> 

p = -[Bsin 2 <j>+CB 2 sin^(()-6 ,2 cos 2 <j))" 5 ]/(7 cos<| 


(2.15) 


(2.16) 


where 


A - 2e I (e ] . + e m ) 

B = e in (B. I + e IH ) t 
C = 2e l£n 
D = 2e. 


E = 2e, 


(2.17) 


The elements of the diagonalized permittivity tensor are 
e, = (1-a) - j'b 
e n = (1/2) if-d) + (J/2) (c-e) 
e ln = [a-(l/2) (c+e)] + 5[^-(l/2)(/+i)] 


where 


a ;■ c !/ 2 w V 

6 " ( 5 “/ /2 “V T5/2 ( “/.V ( 

2 2 

c = [(0 (ti3-iu„) /mV ] C, (m-ui,,) /v 

L p H' m J 3/2 E m 

■ 2 

i = (5» p /2»v m ) C 5/2 (a.-ul B )/v m 


(2.18) 


(2. 19) 


f = (5a) /2wv ) C , (m+to )/v 
J p m 5/2 v E- m 

and v m is the collision frequency associated with the most probable electron 
energy, assuming a Maxwell -Boltzmann distribution of electron energies. The 
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semi-conductor integrals are 


c p (as) - (1/p 


.) 


-u 


2 2 
q 


du 



( 2 . 20 ) 


Although the Sen-Wyller equations provide the necessary tool for 
determining the electron concentration, an explicit solution for N is not 
available because of the great complexity of these equations. In practice, 
electron concentrations are calculated from differential absorption and 
Faraday rotation rates by an iterative computer algorithm [Meehtly et at. 3 
1970] which is based on the Sen-Wyller equations, (2.12) through (2.17) . 
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3. RADIO PROPAGATION AT THE MAGNETIC EQUATOR 
The feasibility of measuring the small differential absorption and 
differential phase rates anticipated at Punta Chilca are discussed in 
Sections 3.1 and 3.2. The rates near the geomagnetic equator are only 0.1 
to 0.01 times those typically observed at midlatitudes. 

3.1 Differential Absorption at the Magnetic Equator 

To ascertain whether or not differential absorption measurements are 
feasible at the geomagnetic equator, radio-wave absorptions anticipated 
at Punta Chilca were calculated by the FORTRAN program ANALYSIS of DA, with- 
out FR, given CF model and the subroutines COEF, FIELD and SENYWL [ Meohtly 
et at. , 1970], A model profile of electron concentration, shown in Figure 
3.1, was constructed from a composite of the profiles of the Croatan series, 
14.228, 14.230 and 14.232 [Meohtly et al. , 1969]. All of these profiles 
are daytime profiles, x (solar zenith angle) 'v 60°. Electron concentra- 
tions at particular altitudes were found by graphical interpolation. Cor- 
responding collision frequencies were obtained from the equation 

V*X:-p C3.i) 

5 - 12-1 

where K equals 7.5 x 10 N ms and p is the neutral atmospheric pressure 
_2 

(Nm ) taken from the COSPAR International Reference Atmosphere Model 
[CIRA, 1965] . For this particular situation, v was calculated at 5 km 
intervals, from 50 to 100 km. Rocket altitudes and velocities were select- 
ed from the trajectory of Nike Apache 14.440. Calculated single mode 
absorptions and differential absorptions are given in Table 3.1. 

Knoebel et al. [1966] estimated that the maximum resolution for dif- 
ferential absorption measurements is approximately ± 0.20 dB s V, or ± 

.133 dB km -1 , assuming a rocket velocity of 1500 ms -1 . Similar estimates 
were reported by Meohtly [1974] , and Ginther and Smith [1975] . A direct 
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TABLE 3.1 

Anticipated single-mode absorptions and differential absorptions. 

Frequency = 2.114 MHz 


ALTITUDE , 
(km) 

O-MODE 

Absorption 

(dB/km) 

X-MODE 

Absorption 

(dB/km) 

DA 

(dB/km) 

65 

.014 

.014 

.0007 

70 

.112 

.128 

.016 

75 

.173 

.220 

.046 

80 

.120 

.171 

.051 

85 

.088 

.128 

. 040 

90 

.149 

.220 

.071 

95 ^ 

.375 

.623 

.248 

100 

Reflection level 
Frequency = 3.145 MHz 


65 

.010 

.010 

.0004 

70 

.065 

.071 

.0059 

75 

.088 

.100 

.012 

80 

.055 

.065 

.010 

85 

.040 

.047 

.007 

90 

.067 

.079 

.012 

95 

.160 

.193 

.033 

100 

.406 

.577 

.171 
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comparison between the estimates of system resolution and the anticipated 
results of Table 3.1 shows that the sensitivity of the experiment is too 
low to measure differential absorption at the geomagnetic equator for solar 
zenith angles of 60° or larger. 

3.2 Differential Phase or Faraday Dotation 

At middle latitudes, the two circularly polarized characteristic waves 
are radiated from the ground at different amplitudes. They combine to form 
an elliptically polarized wave whose axis rotates as the wave propagates 
through the ionosphere because of the changing relative phase of the two 
component characteristic waves. An analogous situation occurs at equatorial 
latitudes. At the geomagnetic equator, the characteristic waves are linearly 
polarized and mutually perpendicular. When these characteristic waves com- 
bine, a resultant wave is produced whose polarization varies with the 
relative phase of the two characteristic waves. Whenever the two character- 
istic waves are in phase or 180° out of phase the resultant wave will be 
linearly polarized. However, when the characteristic waves are out of 
phase the resultant wave is elliptically polarized. The variations, caused 
by the ionosphere, of the polarization of the resultant wave is analogous 
to Faraday rotation at middle latitudes. An illustration of the concept 
of Faraday rotation at the geomagnetic equator is presented in Figure 3.2. 

P. E. Monro [Edwards, 1974] expressed the field vectors for the linear- 
ly polarized ordinary and extraordinary waves at the magnetic equator, and 
derived equations for the rocket receiver output signal, W, as follows: 

= x A sin {(.w-03^/2 )t - k Q Z } - ypA cos{(w-a) ff /2)t - k Q z) (3-2) 

and 

5^ = x. p RA sin {(w+Wd/2) - + y R A cos{ (W+to^/2) t - k^z} (3.3) 


where 


t=l/4 


t=3/8 


t = 0 t=l/8 


t=l/2 


Figure 3.2. Illustration on the concept of Faraday rotation at the 
equator [Rao, 1967]. 
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w = 2 tt/ = propagation frequency (rad s’*) 

oj d = 2tt/ d = frequency differential between the ordinary and extraordin- 

i? ti 

ary wave, often referred to as the reference frequency 
(rad s ’■■*■) 

k 3 k = wave numbers for the ordinary and extraordinary waves (tad m - -*-) 
R = ratio of the ordinary wave amplitude to the extraordinary wave 
amplitude 
A = constant 

x 3 y = unit vectors in the x and y directions, respectively 
p = wave polarization . 

The receiving antenna in the rocket payload is perpendicular to the 
rocket axis. Thus , the unit vector l in the direction of the antenna is 
given by 

l = -x sin in t + y cos w t (3.4) 

s s 

where w = 2 t xf - rocket spin rate (rad s 
s s 

The signal at the output of the antenna is proportional to the scalar 
product of (3.4) with the sum of (3.2) and (3.3). Substituting z = v„t for 
the z coordinate of the antenna into this product yields: 

V = -sin(w-u) /2-fc v )t sin w t - R sin(co+io„/2 -k v n )t sin m t 

H OH S ti X n S 

+ R cos (co+u) n /2 -k vS)t cos oj t - p cos (u)-oi D -7i v D )t cos u> t (3.5) 

ti X ti S ti 0 ti S 

where v ^ is the rocket speed (m/s). Applying the trigonometric identifies 
for the cosine of the sum and difference of two angles to equation (3.5) one 
obtains 

V = cos (oi -(«_/ 2 -7c u+to )t - -hp- cos (oj-a) r) /2-k y -to )£ 

. R . , . ' • 

+ j (1-P) cos (w + u^_/ 2 - k^v t 

•• /"v \ ... ; •' 

+ -■ (l+p) CosCw+w^/2-k^+wp-fr 


(3.6) 
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Equation (3.6) represents the amplitude-modulated signal detected in the 
receiver. The output of the receiver is simply the modulation signal which 
can be extracted from the above equation by converting it to complex ex- 
ponential form and factoring out exp Q/wt) . The modulus of the resultant 
expression is the receiver output W which becomes 

= [jy- (1+P 2 ) + (1-P 2 ) cos 2u ) _t 


W 


+ § (1-P) 2 COS (to^ + ' 2o) s }t v (3.7) 

- j (1+p) 2 cos {u R + (k 0 -k x )v R + 2to s >0' 5 

A second -order binomial expansion and normalization to reduce the dc 

term of this expression to unity yields 

W = 1 + D cos 2to t + E cos {a) E .+a) r ,-2aj }t 
s RFs 


+ F cos {w D +w a +2oj }t - 1/2 DE cos {aj r .-Ho_-4u }t 

zt r S n r S 


(3.8) 


+ ( DE+DF ) cos {co D +a)_}t + DF cos {w n +w ri +4to }t 

n r nr S 


where 


h = = o< 0 -m v r 


D 


E = 


(3.9) 


(1? 2 -1) (1-p 2 ) 

2 (1+i? 2 ) (1+p 2 ) 

fl(l-p) 2 

2 (1+i? 2 ) (1+p 2 ) 

y = -i?(^p ) 2 

2 (1+i? 2 ) (1+p 2 ) 

Since the characteristic waves are linearly polarized at the geo- 
magnetic equator, we set p equal to zero or infinity. 

The principal components of the receiver output signal containing the 

Faraday rotation rate, w n , are 

r 
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E cos <s)^t (3.10) 
F cos (3.11) 
-1/2 DE cos (3.12) 
( DE+DF ) cos w 3 t (3.13) 
DF cos (3.14) 


where 

f l “ V 2 * 5 4 + f F - % 

f 2 = U) 2 /2tt = f R + f F - 2f s 

f 3 = u) 3 /2tt - f ft + f F (3.15) 

'V f 4 m »/**;■? f R + f F + 2f s 

f s = ^ 5 /2tt = f R + f p + 4/ s 

Examination o£ equation (3.8) for the receiver output reveals that the 
components given in equations (3.10) and (3.11) are dominant. Differentiat- 
ing both (3.10) and (3.11) with respect to 5 and setting the derivatives 
equal to zero enables us to obtain conditions on R such that amplitudes of 
terms containing will be extrema. 


dE _ . (l+5 2 )-25 2 _ 

® * 1/2 = 


= 0 


S' 


1-5 


2 2 

2(1+5 ) 


= 0 


1-5 = 0 

5 = 1 


Similarly, 


dF , , n (1+5 ) +25 

3 ? 1/ 4 ‘ , , ^ 2 


5 2 -l 
(1+5 2 ) 2 


= 0 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


(3.21) 
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R 2 - 1 = 0 (3.22) 

R = 1 (3.23) 

By substituting values greater than or less than 1 for R into the 
expression for E and F one can easily verify that R - 1 maximizes E and F, 
That is, the Faraday amplitudes are maximized if equal ordinary and extra- 
ordinary wave amplitudes are transmitted. 

3.3 System Modifications for the Equatorial Radio -Propagation Experiment 

Instrumentation for the propagation experiment at midlatitudes is 
described by Knoebel et at. [1965] . 

With results of Sections 3.1 and 3.2 in mind, the following modifi- 
cations were adopted for the equatorial propagation experiment at Punta 
Chilca, Peru: 

1. Linearly polarized modes of equal amplitudes were transmitted. 

2. The aspect magnetometer sensor and the rocket receiving antenna 
were aligned parallel to each other. 

3. The differential absorption experiment was excluded. 

A system using linearly-polarized antennas for transmission and 
reception, and a frequency differential of 500 Hz between the ordinary and 
extraordinary waves would produce an output at the rocket receiver as shown 
in Figure 3.3. An idealized situation in which the rocket velocity vector 
and the wave normals are in the same direction, and which neglects any non- 
linearities in the rocket -borne receiver is assumed . The axis of the aspect 
magnetometer is aligned with the ferrite rod antenna inside the rocket for 
this example, and in the case of payload hardware, to simplify data inter- 
pretation. Thus, the ordinary wave is sampled during the zero-crossings of 
the magnetometer output, while the extraordinary wave is sampled when the 
magnetometer output is a maximum or minimum as demonstrated in Figure 3.3. 
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In each payload, the aspect magnetometer sensor was aligned with the 2.114 MHz 
receiving antenna. 

Two crystal controlled oscillators, separated by (nominally 500 Hz), 
generate the ordinary and extraordinary wave frequencies. The frequencies 

of the ordinary and extraordinary are given by f - f - 250 Hz and f = f 

o g oo g 

+ 250 Hz, respectively. f Q is the center frequency chosen for a given rocket 
shot. In the case of Nike Apache 14.532, the center frequencies for the 
radio-wave propagation experiments were 3.145 MHz and 2,114 MHz. Each mode 
is transmitted to the antenna dipoles by two 1-kW transmitters. The 
frequency, amplitude, and phase of the voltages applied to the dipoles are 
controlled so that both magnetoionic modes are radiated vertically. A 
block diagram of the system is shown in Figure 3.4. 

The transmitting antenna arrays, shown in Figure 3.5, consist of 
two half-wavelength dipoles arranged at right angles to each other. The 
dipole radiating the ordinary mode was aligned parallel to the magnetic 
field, while the dipole radiating the extraordinary mode was aligned per- 
pendicular to the magnetic field. The dipoles for each array were supported 
by three vertical wooden poles, one pole at each vertex of an isoceles 
right triangle. The poles for the 2.114 MHz antenna array were placed 76 
meters apart and the dipoles were elevated approximately 24 meters above 
the ground. The poles for the 3.145 MHz antenna array were placed 52 meters 
apart and the dipoles were elevated approximately 19 meters above the ground. 
The feed lines of the antenna were connected to coaxial cables originating 
from the transmitters at the ground center of the array. Baiun coils were 
provided for matching the balance-to -ground feed lines to the unbalanced 
coaxial cables. 

The rocket receiving antenna is a magnetic dipole consisting of two 
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Figure 3.4 Block diagram of the modified system. 
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ferrite rods wound with a coil and tuned with a variable capacitor. A one- 
turn shielded symmetrically link couples the antenna to the receiver with 
the proper impedance. A schematic diagram of the rocket receiving antenna 
is shown in Figure 3.6. The antenna assembly is encapsulated in foam and 
mounted within a fiberglass cylindrical section of the payload as shown in 
Figure 3.7. 

The output of the antenna is fed to a transistorized, crystal control- 
led receiver designed and developed by G. W. Henry, Jr. [Edwards, 1975] . 

The receiver uses dual -gate MOSFET devices for the RF amplifier, mixer, and 
three IF amplifier stages. These devices are characterized by a strong 
immunity to overload effects on strong signals, possess a nearly straight- 
line logarithmic AGC control characteristic and are gain-temperature in- 
dependent. The bandpass of the receiver is shown in Figure 3.8 and the AGC 
and detector characteristics versus input signal level are shown in Figure 
3.9. The specifications of the receiver are summarized in Table 3.2. Three 
completed receivers, one with the cover plates removed to reveal the com- 
plete receiver circuitry, are shown in Figure 3.10. The entire receiver 
assembly consists of two printed circuit cards, one containing the RF section 
and one the IF section with input connector, separated by a center metal 
plate. 

3.4 Payload Description 

The payload contained a dc/ Langmuir probe and the radio receivers . 

The Langmuir probe was mounted on the nosecone tip. The receivers and 
antennas for the radio -propagation experiment were enclosed in the fiber- 
glass section of the payloads. The payload also contained the spin magneto- 
meter, a 210 kHz tone range receiver and filter, two quadraloop tone ranging 
antennas and four 30° -turnstile antennas. The four 30° -turnstile antennas 
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Figure 3.9. Detector characteristics versus input Signal level [G, W. Henry, Jr., in 
Edwards, 1975]. 
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TABLE 3.2 


IR-1274 receiver performance data. 


Input : 


Selectivity: 


Frequency 

Level 

Impedance 


-6 dB bandwidth 
-60 dB bandwidth 
Shape factor (60/6) 


Spurious Responses (for 2.225 MHz f ) 
(2.0 to 5.) MHz range) 


2.0 to 5.0 MHz 
-110 to + 10 dBm 
To match ferrite antenna 
(600 ohm nominal) 

2.811 kHz 
4.490 kHz 
1.597 : 1 


2122 kHz 
2156 kHz 
2451 kHz 
2470 kHz 
3135 kHz 


-103 dB 
-103 dB 
-102 dB 
-100 dB 

- 50 dB (Image) 


AGC Characteristics: 

Detector Output Flatness 

AGC Voltage Range: 

AGC Output Impedance 

Detector Characteristics: 

Nominal dc Output 
Linearity 

Temperature Stability 
(-10 to +50°C) 

. Output Impedance 

Power Requirements : 


±0.2 volt, -110 to +10 dBm 
input 

+5.0 volts (minimum signal) 
0.8 volt (maximum signal) 
2.2 kO or kiloohms 


2.0 ± 0.2 volts 
±2 dB, 0.15 to 4.0 volts 
Gain: ± 3 dB 

Offset : ±0.05 volt 

Less than 10 ohms 

+20 to +35 volts @ 55 mA 
-20 to -35 volts @ 15 mA 


Environment : 


Temperature 
Shock and V ibration 


-10 to +50°C 
0 to 15 g, sinusoidal 
30g, random 
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were located at the base of the payload. The two quadraloop tone ranging 
antennas were located directly in front of the turnstile. 

Nike Apache 14.532 included two additional probe experiments for 
measuring electron temperature. These probes were mounted on two deployable 
booms. One of these probes duplicated the boom-mounted probe flown success- 
fully on Nike Apaches 14.475, 14.476, 14.513 and 14.514 [ Schutz et al. , 1975']. 
The other probe was an RF resonance probe supplied by Dr. K. Hirao and Dr. 

K. Oyama from the University of Tokyo. 

The dc/Langmuir probe was modified in this payload to observe the 
fine structure of the electon concentration profile [L. G. Smith in Edwards, 
1975]. The modifications in the instrumentation were suggested by Prakash 
et at. [1972] whose experiment was a direct adaptation of the de-probe ex- 
periment developed by Smith [1967]. In the dc-probe experiment the current 
to the nose tip electrode is measured at a fixed voltage. An electrometer 
with a diode in the feedback loop is used to give a logarithmic scale of 
current over six decades, from lO - ^ A to 10"^ A. The output of the log- 
arithmic electrometer is telemetered to the ground where it is recorded on 
magnetic tape. 

To examine the structure of the ionosphere with greater resolution, 
both in altitude and electron concentration, the ac component of the electro- 
meter output signal was amplified and telemetered to the ground on a 
separate channel. The amplitude of this signal is proportional to the 
fractional change in electron concentration. The circuit of the ac ampli- 
fier used in the Peru launches is shown in Figure 3.11. With this modifi- 
cation we can measure the fractional change in electron concentration to 


0 . 1 %. 



•Figure 3.11. AC amplifier for electron -density fine-structure experiment [L. G. Smith, in 
Edwards, 1975]. 
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4. DATA ANALYSIS 


4.1 An Algorithm for Determining the Spectra of Finite Lengtn Discrete 
Time Sequences 

The problem under consideration here is the determination of the 

frequencies, f. (3.15), of the principal relative maxima, M., of the amplitude 
■ 3 

spectra of finite length discrete time sequences. The discrete Fourier 
transform (DFT) [Bergland, 1969] of a finite length discrete time sequence 
can be calculated by the fast Fourier transform (FFT) [ (International Mathemat- 
ical and Statistical Libraries , Inc. s 1975.] Although most of the properties 
of the continuous Fourier transform (CFT) are retained, several differences 
result from the constraint that the FFT must operate on a limited number of 
discrete samples. The errors most often encountered in using the FFT are 

aliasing and leakage [Bergland, 1969], In order to determine the f.'s 

3 

accurately, aliasing and leakage errors must be minimized. 

The term aliasing refers to the distortion in the amplitude spectrum 
of a signal by the appearance of high frequency components as low frequency 
components. This error occurs as a result of sampling at a rate less than 
the Nyquist sampling rate, f^ . The Nyquist sampling rate is the sampling 
rate at which a continuous function g(t) can be uniquely recovered from its 


sampled version y y ; 

CO 

gm = l gimm-m 

fl* -co 


where T = l/f^ s [Brigham, 1974] . 


git) - T l g{nT) sin 
n=-“> 


In particular, git) is given by 

w Ns (t-m} 

vit-nT) 


C4.1) 


(4.2) 


To eliminate aliasing one simply samples at a rate greater than the Nyquist 


sampling rate (e.g., sample at a rate greater than twice the highest fre- 
quency component present) . 
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The leakage error is related to the way in which a finite-length 
signal is formed from a long-duration signal. Multiplication of the signal 
in the time domain by a rectangular window is equivalent to convolving the 
amplitude spectrum of the rectangular window with the spectrum of the long- 
duration signal. For simplicity, we chose a cosine waveform of constant 
frequency, f , to illustrate the concept of leakage in Figure 4.1. Notice 
that the resultant amplitude spectrum is not restricted to one frequency, 
f , but in fact has a series of small sidelobes. The occurrence of these 
undesirable sidelobes is called leakage. 

Although leakage is an inherent error in the Fourier analysis of any 
finite-length signal, it can be reduced significantly by utilizing a time 
domain window which has small sidelobes in the frequency domain. Several 
examples of windows which have been used by others [Oppenheim and Sohafer, 
1975 ; Blackman and Tukey , 1958] to reduce leakage are the 
Cl) Bartlett window 

L-l 


W{n) = 


2n 

L-l 

2 - 


2n 

ITT 


0 < n < 

— — L 

L-l , t , 
~ 2 ~ — n — 5-1 


(4.3) 


where n and L are integers , 
(2) Hanning window 


A7(n) = y [1-cos (|^)] , 0 < n < L-l 


(4.4) 


(3) Hamming window 


,2irn. 


W{W) = .54 - .46 cos (g^) , 0 < n < L-l 

and (4) Blackman window 


(4.5) 


W{n) = .42 - .50 cos (ll^) + .08 cos (ll2i) , 0 < n < L-l (4.6) 

Jj- 1 — 
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While each of these windows reduces leakage to some extent, a better window 
has been chosen for the present application. 

The Fourier transform of a Gaussian function is another Gaussian func- 
tion. Thus, a Gaussian window eliminates leakage because it has no sidelobes 
in the frequency domain . 

The amplitude spectrum of the rocket receiver signal after windowing 
is Gaussian within a symmetric neighborhood of the / . * s . This is illustrated 
in Figure 4.2. Utilizing this result, E. K. Walton has suggested an al- 
gorithm which accurately determines the frequencies of the principal 
relative maxima of the receiver output containing the Faraday rotation rate. 
Consider the Gaussian curve 

X' = A exp (4.7) 

d 

shown in Figure 4.3, where A and a are constants, and f . is any one of the 

d 

frequencies of equation (3.15). Given three points on this Gaussian curve 


'4-1- *VP. <fk- *'*>> (4... 

where X’^, X' y, X are the amplitudes of the frequency components 

and ’ res P ec tively, and we wish to determine f. . First compute the 
natural logarithm of X'y/X'-^^ and X’^/X'^_^. 


in 


in 


k + 1 


= in 


A exp {-ot( 
A exp 


X' 


k+i 


= -<*( f k -f d r * 


■fc+i ■ g* 


Similarly, 


Aa 


X. 


— * - a <4-4° * ' ,( 4-i 


Dividing equation (4.3) by (4 .4) yields 


Cr = 


- ( 4-4 )2 " oWr**/ 

-t4"4 )2 * ( 4-r4 )2 


(4.8) 


(4.9) 


( 4 . 10 ) 


(4.11) 
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where 


G = 




En c-V r s-_i> 

Expanding the numerator and denominator of equation (4.11) and simplyfing 
yields , 


(4.12) 


G = 


- 2 4.i4 - fj_ « 2 44 
h- 1 - - f k 2 * 2 Vi 


(4.13) 


Multiplying' both sides of equation (4.13) by the denominator on the right 
side of equation (4.13) one obtains 

(4.14) 


g ^-i 2 -4-i ffti* 2 W - h 2 - 2 h>iffh 2 ' 2 Vt 


Solving for f . yields 
0 


„ f k*i 2 - fk - G <4.i 2 - 4 2 > 

4 - - 2 ( 4 - 4 *^ 


(4.15) 


Since the FFT algorithm calculates the Fourier coefficients of the Gaussian 
curve for every integer frequency, equation (4.15) can be simplified. Let 

?'k = k ’ fk-1 = k ~ lr and fk+l = k+1 


then equation (4.15) reduces to 

f - k + k 

T o 2 l+G 

Before the algorithm was applied to rocket data, it was tested on 


(4.16) 


several computer -generated sinusoidal signals of various input frequencies, 

f , and phases, <f> . The specific sinusoids and the results, /., of their 
G 3 

processing are shown in Tables 4.1 and 4.2. The sinusoids listed in Table 
4.2 had been previously processed by an algorithm written by K. L. Miller 
[Edwards, 1975]. Miller's algorithm uses a rectangular window in the time 
domain which does not minimize leakage. Tables 4.1 and 4.2 show that al- 
gorithm (4.14) reduces leakage such that the maximum error in recovering 
the frequency of the input sinusoid is less than ± .0025 Hz. The computer 
program for implementing this algorithm is given in Appendix I. 


I 


TABLE 4.1 

Test #1 o£ algorithm 4.14. 



f * 

f ** 

f.-f 
. 3 J < 

#1 

7.350 

7.35185 

.00185 

#2 

7.350 

7.35162 

.00162 

#3 

7.350 

7.35187 

.00187 

#1 

7.31416 

7.31633 

.00217 

#2 

7.31416 

7.31607 

.00191 

#3 

7.31416 

7.31635 

.00219 


503.5 

503.5 


503.49976 

503.49976 


-.00024 

-.00024 


^Computer generated sinusoids 

Mc N \ 

#1 sin Lhnh" I > N ~ l > 2 , •••» 6000 


(2tf N \ 

#2 sin(^~ f 1 1 ; A 7 = 1, 2, .... 6000 


#3 sin 


2ir/ o ff 


; N = 10, 20, . . . , 6000 


Found by algorithm 4.14 


TABLE 4.2 


Test #2 of algorithm 4.14, 


Phase, 4> 
(Degrees) 

f . * 
3 

f ** 

J 3 

0 

7.00002 

6,978 

10 

7.00002 

6.980 

20 

7.00001 

6.984 

30 

7.00001 

6.989 

40 

7.00000 

6.996 

50 

6.99999 

7.004 

60 

6.99998 

7.011 

70 

6.99997 

7.016 

80 

6.99997 

7.020 

90 

6.99997 

7.020 

100 

6.99997 

7.020 


Computer generated sinusoids 

( 2 *f a N \ 

sin \ 6 ooT 7 + + | N = : ; i* 2 > •• •> 6000 

f = 7.00 Hz 
J o 

* Found by algorithm 4.14 
** Found by Miller algorithm 
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4 . 2 Data Tapes 

Digital data tapes for Nike Apache 14.532 were provided by D. E. 
Hoskinson, of the NASA Wallops Flight Center Telemetry Section. A catalog 
of these post-flight tapes is given in Appendix II. Each tape begins with 
four header records which contain identification information. The next five 
records, of 1005 words each, give the calibrations of telemetry discrimin- 
ators. The subsequent records are data records of 2008 words each. Words 
0006 through 2005 are data samples. Words 2006, 2007, and 2008 express the 
universal time of the first data sample of the next record. The time is 
coded as outlined in Table 4.3. 

The original telemetry data from Peru were recorded on analog tapes. 

The post-flight digital tapes were generated at the Wallops Flight Center 
from the original analog tapes. To maintain accurate synchronization of 
time between the analog and digital tapes, a 100 kHz reference signal, 
recorded on the analog tape during the flight, was used to trigger the 
analog-to-digital converter. By dividing the 100 kHz reference by four, 
each of the five channels could be sampled at a rate of 5 kHz in syn- 
chronization and faster than the Nyquist rate (approximately 1 kHz) 
in a simple way. The circuit designed for this purpose by L. J. Johnson 
of the Aeronomy Laboratory is shown in Figure 4.4. 

4.3 Faraday Rotation 

The Faraday rotation rates were extracted from the radio-propagatiqji 
data by applying algorithm (4.14) to the receiver and reference signals. In 
the processing of these signals a Gaussian window of the form 

Win) = exp(-.56 2 )-Cy (4.17) 

where 

6 = 3(1.0002 - 2!^) ; i = 1 , 2, ..., 2500 


TABLE 4.3 


Coding of time words 


Word 2006 bits 

13 - 16 

9-12 

5-8 

1 - 4 

Wood 2007 bits 

13 - 15 

9-12 
5-7 
1 - 4 

Word 2008 bits 

15 - 16 

11 - 14 

7 - 10 

5 - 6 

1-4 


hundreds of milliseconds 
tens of milliseconds 
ones of milliseconds 
tenths of milliseconds 

tens of minutes 
ones of minutes 
tens of seconds 
ones of seconds 

hundreds of days 
tens of days 
ones of days 
tens of hours 
ones of hours 
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and 


cy = .01111 


was employed to reduce leakage. The factor Cy was introduced to shift the 
Gaussian window downward so that the endpoints would lie on the absccia 
as shown in Figure 4.5. A one-second window was employed to simplify the 
computer processing. 

In Section 3.2 we derived an approximate expression for the receiver 
output under ideal conditions, equation (3.8). The frequencies of the 
principal relative maxima of the amplitude spectrum containing the Faraday 
rotation rate were given by , and Since is independent 

of the rocket spin rate, any change in f will appear only in f 2> 
and jfg. Therefore, provides a reference for identifying fy 

and f- in the spectrum. For example, consider the amplitude spectra of the 
3 MHz receiver output of Nike Apache 14.532 before and after the deployment 
of the boom probes as shown in Figure 4.6. The principal relative maxima 
corresponding to frequencies and f ^ are labeled M^, 

and M ^ in Figure 4.6. Since the deployment of the boom probes decreases the 
rocket spin rate, M and M ^ are closer in frequency to after the 
deployment of the boom probes as illustrated in Figure 4.6. 

The measured quantities jy, f R , and may be combined 

in many different ways to determine fp. For example, 

(4.20) 

(4.21) 

(4.22) 
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Figure 4. 


6 Amplitude spectra of 3 MHz receiver for Nike Apache 
14.532 before and after deployment of boom probes 
(i.e., 51.1 S ET). 
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f i+f 5 

?F 2 ? • 


R 


- f 2+U 

T F 2 


f; 


R 


(4.23) 

(4.24) 


are seven different expressions for calculating f p . Although equation 
(4.18) is the simplest expression for calculating contains a larger 

error produced by cross talk between the data lines. Utilizing either 
equation (4.19), (4.20), (4.21), or (4.22) to calculate requires a 
knowledge of f in addition to / D . Since the receiver and aspect magneto- 
meter outputs contain a large dc component, their amplitude spectra are not 
Gaussian within a 6 Hz symmetric of the spin rate. Therefore, f cannot be 
accurately determined from the lower portion of the receiver spectrum or the 
aspect magnetometer spectrum. While f can be determined by algebraically 
combining or with either or equations (4.23) and (4.24) offer 
a simpler alternative for calculating f„. While equations (4.23) and (4.24) 

n 

are algebraically similar, equation (4.23) contains a large error component. 

This result is illustrated in Figure 4.7. Figure 4.7 is a plot of f versus 

s 

elapsed time, in which f was determined by algebraically combining and 
with and . The large error in illustrated by the large fluctuations 
in f determined by algebraically combining f\. with and is the major 
source of error in equation (4.23). Therefore , equation (4.24) was used to 
calculate f„. 

j F 

The frequencies and the corresponding reference 

frequencies and Faraday rotation rates for the 3 MHz radio-propagation ex- 
periment of Nike Apache 14.532 during the period from 70.5 to 77.0 seconds 
are given in Tables 4.4 and 4.5. No reliable Faraday rotation rates were 
obtained before 72 or after 75 seconds. 

The method of analysis [Meohtly et at . , 1970] based on the complete 


(DE6/S) 



Figure 4.7 Plot of f versus time as calculated j 
v and f $ , s ' 
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TABLE 4.4 



Frequencies f , f 2> 

f A > and f 5 for 

the 3 MHz radio- 



propagation 

experiment 

from Nike Apache 

14.532. 


Elapsed 

time 

(seconds) 

A 


A 

>4 

A 

70.5 

484.1306 

492.41772 

500.41187 

508.66772 

520.77930 

71.0 

484.20313 

492.50122 

500.60815 

508.77734 

520.03027 

71.5 

484.57910 

492.49512 

500.58813 

508.67505 

519.47241 

72.0 

484.35791 

492.49658 

500.52686 

508,74756 

520.57861 

72.5 

484.52246 

492.49536 

500.63306 

508.70752 

519.51074 

73.0 

484.34521 

492.46777 

500.54761 

508.73218 


73.5 

484.35181 

492.47119 

500.51587 

508.77222 

518.97437 

74.0 

484.08081 

492.46265 

500.46265 

508.82715 


74.5 

483.00000 

492.45947 

500.58203 

508.81689 


75.0 

485.54980 

492.45850 

500.47632 

508.77002 


75.5 


492.27661 

500.49585 

508.72559 


76.0 

485.90576 

491.86963 

500.59277 

508.72144 

520.98535 

76.5 


492.53223 

500.56079 

508.82617 


77.0 


492.19604 

500.47046 

508.80420 


77.5 


492.41162 

500.38013 



78.0 


492.39453 

500.69238 

508.90430 
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TABLE 4.5 


Reference frequencies and Faraday rotation rates for 3 MHz 
radio -propagation experiment from Nike Apache 14,532. 


Elapsed 

Time 

f R (Hz) 

f 2 4 ’^ 4 

2 - h M 

fp (Hz) 

fp (deg/ 
sec) 

70.5 

500.51147 

500.54272 

.03125 

11.2500 

71.0 

500.57861 

500.63928 

.06067 

21.8412 

71 . 5 

500.59106 

500.58509 

-.32551 

-117.1836 

72.0 

500.58789 

500.62207 

.03418 

12.3048 

72.5 

500.58618 

500.60144 

.01526 

5.4936 

73.0 

500.58569 

500.59998 

.01429 

5.1444 

73.5 

500.58472 

500.62171 

.03699 

13.3164 

74.0 

500.58105 

500.64490 

.06385 

22.9860 

74.5 

500.57935 

500.63818 

.05883 

21.1788 

75.0 

500.57739 

500.61426 

.03687 

13.2732 

75.5 

500.57422 

500.50110 

-.06322 

-22.7592 

76.0 

500.57056 

500.29554 

-.27502 

-99.0072 

76.5 

500.56665 

500.67920 

.11255 

40.5180 

77.0 

500.56470 

500.50012 

-.06458 

7 23.2488 
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Sen-Wyller equations without approximations [Sen and Wyller, 1960] was 
employed to analyze the Faraday rotation data from Nike Apache 14.532. The 
procedure for analyzing is outlined in Figure 4.8. An initial value for 

the electron concentration is assumed. From this initial value of electron 
concentration, the values of the refractive indices for both magnetoionic 
modes are calculated. The refractive indices are used to obtain a calcu- 
lated value of f corresponding to the initial value of electron concentra- 
tion. The electron concentration is repeatedly corrected until the calculat- 
ed value of f matches the experimental value. The calculated value of f 
t . . • £f . 

and the experimental value of f usually agree to within ± .1% after five 

b 

or six iterations. 

To solve for the electron concentration when only Faraday rotation 

data are available, a model for the collision frequency must be used (i.e., 

equation (3.1)). The pressures were taken from the COSPAR International 

5 -12-1 

Reference Atmosphere Model [CIBA, 1972] and K was equal to 7.5 x 10 N m s 
The electron concentrations obtained from the analysis of the 3 MHz radio 
data of Nike Apache 14.532 are given in Table 4.6 and are plotted in Figure 
4.9. The error bars correspond to an uncertainty of ± 3 deg s ^ in the 
Faraday rotation rate obtained by applying equation (4.22) to the data 
between 35 and 52 km under the assumption of free space conditions. 

4.4 Probe Current Calibration 

Variations of electron concentration with altitude are believed to be 
accurately represented by changes in probe current [MechtZy et at. , 1967] . 
However, exact equations giving electron concentration as a function of 
probe current, collision frequency, rocket velocity, and other parameters do 
not exist. For this reason, dc -probe current, I, is calibrated by values 
of electron concentration, N, from the radio-propagation experiment. The 
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TABLE 4.6 

Electron concentration calculated from 3 MHz 
radio data of Nike Apache 14.532 


Center 

Time (s) ff(m - 3) 

72.0 1.484277 x 10 10 

72.5 7.311761 x 10^ 

73.0 6.887715 x 10 9 

73.5 1.596224 x 10 10 

74.0 2.457148 x 10 10 

74.5 2.314058 x 10 10 

1.605942 x 10 10 


75.0 
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N/I ratio is calculated for each value of N from the radio data analysis 
and from corresponding weighted mean values of I. The values for N, I, and 
N/I for Nike Apache 14.532 are tabulated in Table 4.7. These N/I values are 
plotted as a function of altitude in Figure 4.10. 

The altitudes of the 2 MHz and 3 MHz ordinary mode reflection points 
are derived from the times at which the AGC voltages of the rocket receivers 
indicated a loss of signal. 

The reflection of the 2 MHz ordinary mode occurred at 77.85 S ET which 

corresponds to an altitude of 100.6 km. The critical value of N is 5.54 x 

4 -3 

10 cm at this point. 

The reflection of the 3 MHz ordinary mode occurred at 111.57 S ET which 
corresponds to an altitude of 139.1 km. The critical value of N is 1.23 x 

5 -3 ’ 

10 cm at this point . 

Past experience shows that N/I is nearly constant between about 100 

and 150 km. Therefore , the calibration between .105 and 120 km is set at the 

9 -3 

value determined by the 3 MHz ordinary mode reflection point, 5.6 x 10 cm 



Past experience shows that a transition in N/I occurs’ at the steep 
gradient in electron concentration at the mesopause. The point (2.8 x' 

9 -3 i 

10 cm A , 86/5 km) corresponds to the intersection of a plateau and 
a steep gradient of the probe current. Between 86.5 and 105 km, the equation 
for the straight line passing through the 2 -MHz ordinary mode reflection 

9 -3 -i 

point and the point (2.8 x 10 cm A , 86.5 km) is 

tf/JOs) = 2.8 x 10 9 exp , 86.5 <_ s < 105 km. (4.25) 

The error bars in this region are extrapolations of the uncertainty of ± 3 
deg s" 1 previously determined for the Faraday rotation rates under the 


TABLE 4.7 


Altitude N 

(km) (m _ 5) 


92 .795 

1 .484277 

X 

ioio 

93.473 

7.311761 

X 

10 9 

94. 149 

6.887715 

X 

10 9 

94.822 

1.596224 

X 

10 10 

95.493 

2.457148 

X 

10 10 

96.162 

2.314058 

X 

10 1Q 

96.828 

1.605942 

X 

10 10 


N/I 

(m“^A - -*-) 
7.9586 X 10 15 
2.9626 x 10 15 
2.0999 x 10 15 
3.7070 x 10 15 
4.5085 x 10 15 
.6526 x 10- 5 3.5459x1015 

.7474 x 10-5 2.1487 x 10 15 


The N, I and N/I values from Nike Apache 14.532. 

Weighted I 
(A) 

.1865 x 10“ 5 
,2468 x 10-5 
.3280 x 10" 5 
,4360 x 10-5 
.5450 x 10~ 5 
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assumption o£ £ree-space conditions between 35 and 52 km. 

9-3-1 

The constant value 2.8 x 10 cm A is assumed as the calibration 
below 86.5 km. 

For each second of the flight, ten values of electron probe current 
are multiplied by the appropriate values of N/I from the smoothed curve of 
Figure 4.10. Since the present method of data processing is based on sampl- 
ing of the probe current at equal intervals in time, the points are unevenly 
spaced in altitude because of the changing velocity of the rocket. K. L. 
Miller's program [Edwards, 1973] was used to interpolate electron concentra- 
tions at 0.1 km intervals. The final results are presented in Chapter 5. 
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5. CONCLUSIONS 

The final profile of electron concentration from Nike Apache 14,532 is 
shown in Figure 5.1. Corresponding numerical values of electron concentra- 
tion are tabulated in Table 5.1. 

The model profile, used to calculate the anticipated differential 
absorption rates in Section 3.1, is also shown in Figure 5.1. 

The profile labeled K14 in Figure 5,1 is a profile obtained by Kane 
[1974] at a solar zenith angle of 14 degrees from rocket measurements of 
radio-wave absorption conducted : from the equatorial site at Thumba, India on 
March 8, 1968. Both rockets, 14.532 and K14, were launched on one of the 
five quietest days of their respective months, in terms of magnetic activity 
and disturbances. Each day was designated "QQ" by World Data Center A for 
Solar-Terrestrial Physics- [ Lincoln , 1968; 1975]. 

A comparison of these profiles reveals that: 

1. There is good agreement between the model profile and the profile 
of Nike Apache 14.532. This observation supports the conclusion 
that for a solar zenith angle of 60°, differential absorption at 
the geomagnetic equator is too small to measure by present tech- 
niques. (cf. Section 3.1). 

2. All three profiles exhibit a plateau of electron concentration. 

The plateau in electron concentration occurs at a higher altitude 
in the model and 14.532 profiles than in K14. 

3. The electron concentrations measured by Kane are larger by a factor 
'X of approximately 4. ■ 

4. A steep gradient in electron concentration, beginning at 86.5 km, 
is evident in the profile of 14.532. 

5. The ratio iy?CK14)/iy)(14.532) = 3.4, where iy? is the value of 



10° I0 1 I0 2 I0 3 I0 4 I0 5 I0 6 


ELECTRON CONCENTRATION (CM' 3 ) 

Figure 5.1 Profiles of electron concentration. The inflection points, denoted 
IPT, for 14.532 and K14 are (3.23 x 10 2 cm- 3, 86 km) and (1.10 x 
lO^ cm“3, 82.7 km), respectively. 
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the maximum ionization in the D region, just below the steep 
gradient, 

6. No sporadic-E layers are observed. 

, A characteristic feature of daytime profiles of D-region electron con- 
centration, N, is a gradient just above an inflection of N at altitudes 
ranging from about 80 to 86 km [Mechtly and Bilitza , 1974] . Mechtly and 
Bilitza showed that N^D , the value of N at the inflection point, is pro- 
portional to the square of the cosine of the solar zenith angle, y 

N D. « (cos x) 2 (5.1 

when the D region is unperturbed. 

By equation (5,1), the ratio of N^D for a solar zenith angle of 14 
degrees (i.e., K14) to N D for a solar zenith angle of 60 degrees (i. e. , 
14.532) is: 


N m D( K14) 
i\M?(14.532) 


(cos 14°)' 
(cos 60°)' 


= 3.8 


This compares with 3.4 calculated from the values of N m D shown in Figure 


Thus, profiles K14 and 14.532 support the conclusion that ff Z>«(cos x)' 


when the D region is unperturbed. 
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C 

C 

C 


979 

978 


C 


0 

43 


:i:WK MUST BE FOUND I N CORE ON A BOUNDRY DiUibiBLb B' 6 
INTEGER *4 I WK < 18050 > WILL ONLY WORK HALF THE. lAMt • . 

THIS IS BECAUSE IT IS USED AS A STORAGE AREA FOR A REAL 

$4 VAR I ABLE AN H- i 

REALMS' IWK ( 9025) 

INTEGERS ARRAY 032:1.28) 

REALMS DATA (6000) 

COMPLEX* :l. 6 GAMN 

REAI...Y 4 DMAS ( 3000 ) v FREQ ( 200 ) 

INTEGER# 4 DAT AID! 20 ) 

EQUIVALENCE ( DATA yDMAG) 

NAMEL I ST/PR AM8/NRG0 y NSNGO y NELT y NELN 
ASSIGN 999 TO KILL 
READ 979) DATA ID 
FORMAT (5Xy204//4) 

WRITE (6 ,,578) DATA ID 
F 0 R M A T C 5 X y 2 0 A 4 / / ) 

CALL TP OP I Z ( 12) 

NR GO IS FIRST RECORD NUMBER 

NSNGO IS FIRST SCAN NUMBER IN RECORD 

NELT IS ELEMENT NUMBER IN SCAN 

NFI N IS NUMBER OF ELEMENTS TO BE US'l D * ... I 

WARNING •••• USER MUST COMPUTEALL OF THESE VALUES GAKEELLLr 

PLOTTING SEQUENCE . ********************* 

MH Pi. 30 .. ... 

GENERATE FREQ* ARRAY (MINUS TO MOVE DOWN 1 HE PAbE 

D I MENS I ON SC A (2 ) y SCF ( 2 ) 

BO 8 NF=“.I. y ISO 
FREQ ( NF ) ~ :I. -NF 

FIRST DRAW FREQ* SCALE ... Tjimi ' , VTr ,, 

UTART AT ZEROv THEN GO AT 20* PER INCH(EUR b INCH H/%.1. o>; 

SCF ( 1 ) “"0 * 0 
SCF (2) ™20. 

CALL PLOT ( * 5 y * 5 y "-3 ) . „ 

G ALL CCP5AX ( 0 ♦ y 9 V y / LOW FREQ ( HZ ) 7 » ••••13 y 3 * y "•90:* yGCP i ; 

CALL CCP5AX(0* rd* y 'FREQUENCY (HZ) ' y ••-14 y 6 * y ~‘V0 * > SUE ) 

READY SCALE FACTOR FOR FREQ ARRAY 

SCF ( 1 ) :::: ••••ISO.O aV .:< ' t 

SET UP TIME AXIS* (NCRV « # OF PLOTS 10 MAKE) 

\ SCA ( I ) "" 1 * 

SCA ( 2 ) ™:t. ♦ T ; 7;TV_,j/<\ 

CALL CCP5AX ( 0 * yO* y 'TIME (RELATIVE) 7 y -<.4 y D^ST y 0 . ,SCA> 

CAI I GCP5AX( 0 * y 9 ♦ y 7 TIMET RELATIVE ) 7 y 14 y D.I.S (yO.y&CA) 
PLOTTING SEQUENCE YYYYYYYYYYYYYYYYYYYY 
CONTINUE -c .; v V 

READ ( 5 v PRAMS ) GRjQjjy* T & A 

WRITE ( 6 y 889 ) NR GO y NSNGO y NELT y NELN OP p 0 Op „ IS 

lt DUALITY 


i 


I 


I 


I 


I 
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889 

n 

w 


40 

41 


:l. 0 4 
C 



10 

C 



10 


i::; 

\.j 


FORMAT ( IX » 91 10) 

FORMAT DATA CARD XPRAMS NRGQ=4000 » NSNGQ«20 1 ETC y SEND 

IF<NR0Q*EGUO>. 00 TO 42 
T A PE H A N D L I N G R 0 U T I N E S 
I F ( N R 0 Q , E Q * 1 ) C A L L T P B 8 R Z ( 1 2 > 

:i:f<nrgo*eq*o) go to 42 
iftnrgo.eq.i ) go to 41 

KSKXP - NR GO -l 
DO 40 N-I y KSK XP 
CALL TPFSRZ ( 12 ) 

CONTINUE 

FNTOPeFLOAT < NELNTNSNGO > /400 . 

N T 0 P - X F I X < F N T 0 P ) T 1 

KNBNGD-I F XX ( 400 ♦ # < FNT0P1FL0AT ( l-NTOP ) ) ) 

DO 50 N-lyNTOP 

M - C N •••• 1 ) $ 2 0 08 T 1 

CALL TP GETZ < X2> ARRAY ( M ) ) 

CALL TPGHKZ ( 12 y NT APE y KILL ) 

IG0-MT1999 
XSTOP » X GOT 8 

WRITE < * 104 ) < ARRAY < X ) y X-=XGO r XSTOP ) 

FORMAT (IX s- 9Z10) 

NT A PE - # OF BYTES TRANSFERRED 
CONTINUE 

WRITE ( 6 9 102 ) ARRAY (2006 ) y ARRAY C 2007 ) , ARRAY (2008 ) y NT APE 
FORMAT ( 1H0 i ' TIME CODE Ts-TZlO) 

MOO 5-2005 

NA*NSNG0#5TNELT ' 

D O 10 N-l y NELN 
DATACN)- ARRAY (NA) 

N A -NAT 5 

XFCNA.LE. NO 05 ) GO TO 10 
N A -NATS 

NOO 5 - N 00 54 2 0 08 ; 

CONTINUE 

AT THIS POINT s’ “DATA" IS READY FOR WINDOWING 
NELN2-NELN/2 * 

FLN2-NELN2 

STT=1»T1, /FLOAT (NELN) : 

GBASE - 0.01 111 ( 

DO 20 K-l s’ NELN 2 
DEL :-B.:3*6*< ; STt.-PLGAtXK>:/FLN2:)- 
6-EXP (•••() . 5HfDEL>l<DEL) “"GBASE 
DATA(K) - G)KDATA(K) 

KTP-NELNHLTT T. ; . 

DATA (KTP) - GYDATA ( KTP ) 

CONTINUE 

NOT NEEDED NOW REMOUE CODE COMPLETELY LATER 

IF (NORM vNE » 0) GO TO 7 

XG0-NELN2 - 100 

XSTOP - X (3 0 T 199 

WRITE ( 6 y 105 ) ( DATA ( X ) y X-XGOy XSTOP) 

FORMAT ( 1 0 ( 1 X y E9 . 2 ) ) - 

CONTINUE . -tF-/' ■ '■ / T " 


C A L L F F T R < D A TA y G A M N v N E L N y i U K > 

C AT THIS POINT ? "DATA “ CONTAINS THE COMPLEX SPECTRUM* 

C WRITE IT OUT 

C CONVERT TO MAGNUTUDES. ************** 

DO 30 N-l > NELN2 
NT » m2 

DMAS ( N ) =DSQRT ( DATA ( NT ) *DATA ( NT > +DAT A < NT-1 ) :*DATA ( NT- 1 ) ) 

30 CONTINUE 

C WRITE IT OUT IN SELECTED CROUPS #*######## 

DO 31 N-l y 4()y 10 
NFREQ ■* N-l 
NTP « N + 9 

3 1 WR ITE ( 6 t 1 03 > NFREQ y ( DM AG ( K > y K~N y NTP ) 

103 FORMAT ( IX yI9y 10 (IX y E9 *2) ) 

DO 32 N - 451 y 551 y 10 
NFREQ-- N-l 
NTP a N •}• 9 

32 WRITE (6-y 103 ) NFREQ y (DMAG( K ) yk'-N y NTP ) 

C ■' , ; 

CJK3k*#}|c^FIND FREQUENCY OF PEAK#*###* 

Kl-6 
K2-T6 
MAX- 5 


33 AMAX-0 ♦ 


DO 52 K-K'l y K 2 

IF ( DMAS < K ) * LT * A MAX > GO TO 52 
MAX-K 


AMAX-DMAGCK) 

52 CONTINUE 

CCCCCC COMPUTE FREQ OF SIG * ALGORTHIM «« 
MAXM-MAX-1 
MAXP-MAX+i 


106 




RR-ALOG ( AMAX/DMAG ( MAXP ) ) /ALOG ( AMAX/DMAG ( MAXM ) ) 
FMAX1-MAX- 1 

FSTG-FMAXF. 5# < 1 *0~RR >/l ♦ OfRR ) 

WRITE (6 y 106>FSIG 

FORMAT ( /y ' FREQ* OF PEAK - ' »F10.5> 

K 1-4 80 

XKM1^DMAG(476) 

XKM “DM AG (477 ) 

XK-BMAG(478) 

XKP-DM AG ( 479 ) 


DO 54 K-“ K1 y K2 

IF ( XK ♦ GT ♦ XKM .AND. XKM . GT . XKM1 .AND ♦ XK . GT . XKP ) GO TO 

XKM1-XKM - 

XKM-XK 

XIO-XKP 

XKP-DM AG ( K ) 


GO TO 54 
MAXM-K-3 

RR-ALOG ( XK/XKF' ) /ALOG ( XK/XKM) 
FSIG-MAXMf .5* ( 1 ♦ O-RR ) / < 1 . O+RR ) 


PAqP to 

OP POOR ' lb 
uu QhALiTv 




] 


I 


1 


WRI TE ( 6 f .1.06 ) FSIG 
(30 TO 56 
CONTINUE 

: DANGER<<<< DATA CANT!) GIVE DIVIDE BY 0 
OF FREQ & COMPUTATION ON SPECTRAL PEAKS 

P L OT T I N (3 S E Q U EN C E ******************** 

SET UP SCALING FOR ALL PLOTS 
SCAT 1 ) :s: l * 

I F ( MELT ♦ EQ . 5 ) SC A ( 1 ) “0 ♦ 0 
SC A ( 2 ) «2 * 

CALL XI..OGZ C DM AG » FREQ > 60 , If SCA , SCF ) 

WR I TE ( 6 y 1 0 1 ) SCA < 1 ) * SCA ( 2 ) 

FORMAT < IX *2E 15*5) 

C A L L X L 0 G Z < DMAG < 441 ) v FREQ < 60 ) * 120 * 1 * SCA , SCF ) 

CALL PLOT Cl ♦ > 0 ♦ f-3 ) 

READY FOR NEXT XLOGZ CALL. ♦ 

PLOTTING SEQUENCE ******************** 

READ NEXT DATA CARD IF SCAN # » ~ Z * D° OTHER NELT 

♦ ♦♦IF 8CAN#*0» READ NEXT NEL.N ELEMENTS* 

GO TO 43 

WRITE ( 6 , 998 ) NT APE 

FORMAT C ' TAPE ERROR 7 *110) 

STOP 

END 




APPENDIX II 

Nike Apache 14.532 Post-Flight Data Tapes 

Tape #3, Receiver #1 

File 1 Probe Log 

File 2 Extraordinary Power #1 

File 3 Extraordinary Power #2 

File 4 Ordinary Power #1 

File 5 Ordinary Power #2 

Tape #4, Receiver #1 

File 1 Receiver Modulation #1 

File 2 500 Hz Reference #1 

File 3 Receiver Modulation #2 

File 4 500 Hz Reference #2 

File 5 Magnetometer 

Tape #11 j Receiver #3 

File 1 Probe Log 

File 2 Probe Fine Structure 

File 3 Receiver AGC #1 

File 4 Receiver AGC #2 

Tape #12, Receiver #1 

File 1 RF Probe Signal 

File 2 Boom Probe Log 

File 3 Boom Probe Linear 

File 4 RF Probe Monitor 



File 5 


Magnetometer 


